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Abstract 

We identify conditions for the presence of negative specific heat in non- 
relativistic self-gravitating systems and similar systems of attracting par- 
ticles. 

The method used, is to analyse the Virial theorem and two soluble 
models of systems of attracting particles, and to map the sign of the 
specific heat for different combinations of the number of spatial dimensions 
of the system, D(> 2), and the exponent, 0), in the force potential, 
<f> = Cr". Negative specific heat in such systems is found to be present 
exactly for v = — 1, at least for D > 3. For many combinations of D and 
v representing long-range forces, the specific heat is positive or zero, for 
both models and the Virial theorem. Hence negative specific heat is not 
caused by long-range forces as such. We also find that negative specific 
heat appears when v is negative, and there is no singular point in a certain 
density distribution. A possible mechanism behind this is suggested. 

PACS numbers: 04.40; 05.90.+m 
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1 Introduction 

In the area of self-gravitating systems, the literature often gives the impression 
that the long-range properties of forces are the cause of negative specific heat in 
such systems. Here, we will show that only some systems affected by long-range 
forces exhibit negative specific heat. 

We limit our investigations to potentials of the form 



0(r) = Cr v (1) 

and of the form 

<f>(r) = d\nC 2 r (2) 

where r is the radial coordinate in a space with D dimensions (D > 2), v an 
integer (y ^ 0), and C, C\ and C2 constants (only depending on D and v). 
Observe that gravitation is described by 0) for D ^ 2, and by 10) for D = 2. 
We will, however, to extend the analysis, regard both 0) and (J5J for any D > 2. 
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We can use spaces with D > 2 and v ^ as our theoretical "test bench" 
to investigate the cause of negative specific heat. There is, for instance, no a 
priori reason why long-range forces should result in negative specific heat in 
three dimensions, but not in other number of dimensions. If this happens only 
in three dimensions, there is probably another mechanism laying behind, that 
is more relevant to explain negative specific heat. 

Before investigating negative specific heat, we will define the concept of long- 
range and short-range forces. Let us regard a continuous medium with constant 
density, and investigate from which areas of the medium the main part of the 
potential energy of a particle embedded in the medium comes. Short-range 
forces are, for potentials of the form 0), then characterised by 

D + v < (3) 

and long-range forces by 

D + v > (4) 

In the case where D + v = 0, the dependence of distance for the total potential 
energy of the particle is logarithmic, and the force is then both short- and long- 
range. For potentials of the form J5J), the corresponding criteria implies that we 
have a force that is both short- and long-range. 

The definition of specific heat at constant volume is 

cv . m^w <5) 

where "<>" expresses a time average. Negative specific heat was first inves- 
tigated by Lynden-Bell and Wood and Thirring [2J. It is also described by 
Hut [Hj, Lynden-Bell @], Lynden-Bell and Lynden-Bell [H] and Padmanabhan 0. 

The Virial theorem in its simplest form, without external pressure, applied 
to a system of particles interacting with potentials of the form reads 

< E >= ^ < K > (6) 

Using (0 and the assumption a <^> |y > o, which is valid for many systems, 
we obtain that negative specific heat appears exactly for v = — 1, independent 
of the number of dimensions, D(> 1), of the system. In three dimensions 
v > — 2 corresponds to long-range forces. So, for three-dimensional systems 
(and negative v), negative specific heat seems to be correlated with the long- 
range nature of the forces. The Padmanabhan model and the Lynden-Bell model 
(see below) were invented to verify that v = — 1 gives negative specific heat in 
three-dimensional systems. When considering systems with D > 3, it becomes 
apparent that the domain (in D and v) of negative specific heat predicted by 
the Virial theorem, is just a small part of the domain with long-range forces. We 
will show that generalised versions of the Padmanabhan model and the Lynden- 
Bell model (both with D > 2) give the same domain of negative specific heat 
as the Virial theorem, for D > 3 and any 0). For systems with particles 
interacting with potentials of the form J5J), and without external pressure, the 
Virial theorem reads 

< K >= constant (7) 

That is, energy supplied only contributes to the potential energy of the system. 
Here we clearly have a <^> | y = 0, so this is one of those systems where the 
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spatial volume is important in the contribution to the temperature. We will 
show that, for this kind of system, the two generalized models give Cy > for 
all D(> 2). 

The result of the analyses of systems with potentials of the form (JJJ is 
mapped in fig. 1. The border between long- and short-range forces is marked. 
Gravitation obeys, for D ^ 2, v = 2 — D, and is therefore, according to (@J, 
a long-range force for D ^ 2. Only for v = — 1 we have significant negative 
specific heat in our models. Only the Padmanabhan model has a non- negative 
value on Cy for v = — 1, namely when D = 2. 

It is, however, as often suggested, reasonable to believe that long-range forces 
in a system is the cause of non-extensivity of the system. With non-extensivity, 
we mean that if two originally separated subsystems with energy E\ and Ei 
are combined, their total energy will not be E\ + E%. For systems defined by 
short-range forces, the interaction energy between subsystems becomes negligi- 
ble for large enough subsystems, since the interaction energy scales as the area 
of the subsystems, and then we have an extensive system. For subsystems in- 
teracting by long-range forces, the interaction energy remains significant, since 
the interaction energy scales as the volume of the subsystems, and we then have 
a non-extensive system. Oppenheim 7 investigates the statistical physics of 
systems with long-range interactions. 



2 Two soluble models 

The models under study are the Padmanabhan model and the Lyndcn-Bell 
model. We are interested in the interval of E where the system may have 
negative specific heat, that is, where no inner or outer cutoff affects the system. 
When regarding the outer cutoff, this implies E < for systems with potentials 
of the form (JIJ, and any E for systems with potentials of the form J3J). (We chose 
the radius of the cutoffs so that they affect the system as little as possible.) In 
their original version, these models involve an inner cutoff, with radius a, that 
keeps the phase-space finite. This results in a region in E with positive Cy, 
even for systems where there is also a region where there is negative Cy. Here, 
we want to let a — > 0. This admits a study of the effect of the potentials 
under consideration, without (model dependent) influences from this cutoff. As 
we will see, it also results in constant Cy in the interval of E for which the 
solutions are valid. Therefore, results for Cy that in the original models were 
represented by functions of E, can here be represented by a single number. The 
limit a — > can also be regarded as an idealization of the case where there is 
an inner cutoff with a small radius, but a radius so close to zero, that the value 
on Cy does not significantly differ from the one obtained in the a — > limit. 
To obtain an unambiguous value on Cy when a — > 0, we have to specify that 
other parameters do not go to some limit when a does. This especially applies 
to E for which we define \E\ < Eq, and to N (the number of particles in the 
Lynden-Bell model) for which we define \N\ < Nq, where Eq and No are some, 
eventually very big, numbers. 

The Padmanabhan model does not give well defined values on Cy for D = 1, 
and the Lynden-Bell model is not meaningful for D = 1. Therefore, we limit 
our investigations to D > 2. 
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2.1 The Padmanabhan model 



This model is presented in Padmanabhan [jy. It contains two particles, each with 
mass to, that attract each other. They are contained in a spherical container 
to obtain ergodicity. In this essay, we generalise the model to D dimensions 
(D > 2) and a potential of the form with arbitrary v (y 5^ 0), and where C 
and v have the same sign. The Hamiltonian is 

H(f , q ,Pi,qi) = ^_ + _EL + mCr » (8) 
2too 2toi 

where Too is the two particle's total mass (= 2m), mi their reduced mass (= 
Po the conjugate linear momentum of their centre of mass, pi their relative 
conjugate linear momentum, and r the distance between the two particles (= qi). 
qo is the position of the centre of mass. The volume of the phase-space is 

(0) (E) = C f dr r D - x (E - mCr u ) n ° (9) 
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where no = D — 1, and r max is the limit on r set by the total energy of the 
system, E. The number a represents an inner cutoff that keeps the phase-space 
volume, g(°\E), finite. 

We calculate the specific heat according to 

9 {1) {E? 

° V = k g^(E) 2 -g(o)(E)g^(E) (10) 

where g^(E) and g^(E) is the first and second derivatives of g^(E) respec- 
tively. From now on, we restrict the calculation to the case D > 3. The case 
D = 2, we consider later. The formulas for g^ 1 ' and g* 2 ) differ from g^ with 
other exponents (ni = D — 2 and 712 = D — 3) and with other constants cor- 
responding to C . It follows from the properties of (|10|l . that multiplication of 
gV) with a factor ^ 2+73 ™ 3 arbitrary, but not dependent of j) do not change 
Cv- We let a-tO. We can use this to simplify (fTTT|) . 

The integral in can be solved by partial integration. For D + iv 7^ 
(i = 1, 2, . . . , -D — 3), the result looks like 

&rr D - x (E - mCr») n ' = [r D (E - mCr v ) n i} r a """°+ 



D(D + v) 



[r u+v {E -mCr v ) r 



% -( % --1Kto 2 C 2 g+2 _ mCr u )nj -2 ]r 
D(D + v){D + 2is) 1 y ' ia 

, n A n j - 1) ■ ■ ■ ■ • 1 • v^m^C^j r n^» ]f , 
D(D + v){D + 2v) ■ . . . ■ {D + njvy Ja 



(11) 



When a — > some terms will completely dominate in this expression. When D + 
iv = for some i (i = 1, 2, . . . , D — 3), terms with logarithms of r appear, some 
of them dominating. When only considering dominating terms, and multiplying 
the gU>:s with allowed factors, we get simpler expressions for Cy. 
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First, we consider D + rijV > for all j. We call this case A. We have 
(o) I 

(i) ^ I (\2\ 

y D+(D-2)v \ 1 ^> 

„ 1 

where the "~" sign stands for equality after multiplication with an allowed 
factor and after taking the limit of a. For the specific heat, we have 

,D+(Z>-1> 

C v = k i — D + n jV > (13) 

v 

This expression is positive for positive v, and for negative v negative in the 
complete interval where it is valid, namely for v = — 1. 

We then consider D + rijU < for all j. We call this case B. We have 

„(0) ^ J_ i Cgrj> + {D-l)(D-2)u 2 (D-l)(D-2)-..,l-u D -\ 

y D ' D(D+v) ' D(D+u)(D+2v) Z5(L>+i/)-...-(_D+(D-l)i/) 

»W ~ m IV 1 \ ( D - 2 > I (g-2)(g-3);/ 2 (J-2)(g-3)-...-l^ 13 " 2 \ 

{/ l- 1 ^ J-Ad "I" D(D+u) D(L>+^)(D+2i/) D(_D+i/)-...-(D+(_D-2)i/) 

fl( 2 ) ~ I'D 9Vn IV 1 I ( D - 3 > | (J'-3)(g-4),y 2 
g y ' ~ - Z){D - L){p + D ( D+v) + D{D+v){D+2 u) +■■■ 

(D-3)(D-4):.,l-v D - 3 x 
_D(D+^)-...-(_D+(D-3)^) / 

(14) 

Then, we will consider the case where there is one i (i = 1, 2, . . . , D — 3) for 
which D + iv = 0. We call this case C. We obtain 

fl <°> ~ 1 

ff (i)^£)_j_l (15) 
</ 2 ) - (D-i-l)(D-i-2) 

For the specific heat, we have 

C v = k(D-i-l) 3i:D + iv = i = 1, 2, . . . , D - 3 (16) 

This expression is positive in all intervals where it is valid. 

Now, we investigate the cases for which none of the criteria for case A, B or 
C is fulfilled. We call this case D. We can conclude which g^'-.s that dominate 
in ljl()|) when calculating the specific heat. Depending on which g^'-.s that 
dominate, the specific heat becomes either or k. 

The case D = 2, which we excluded earlier, gives Cy = k for all 0). 

In the table below, numerical values on Cy in terms of units of Boltzmann's 
constant, k, are presented. The type of case is also given. 
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V 




D 






2 


3 


4 


5 


6 


7 


8 


9 


3 




1 


3.0 (A) 


4.3 (A) 


5.7 (A) 


7.0 (A) 


8.3 (A) 


9.7 (A) 


11 (A) 


2 




1 


O.O (A J 


O.U (AJ 


0.0 (AJ 




y.o [A.) 


11 ( A ) 


12 (A) 


1 




1 


5.0 (A) 


7.0 (A) 


9.0 (A) 


11 (A) 


13 (A) 


15 (A) 


17 (A) 







X 


X 


X 


X 


X 


X 


X 


X 


1 




1 


-1 (A) 


-1 (A) 


-1 (A) 


-1 (A) 


-1 (A) 


-1 (A) 


-1 (A) 


-2 




1 


0(D) 


1 (D) 


1 (D) 


2(C) 


2.5 (B) 


3(C) 


3.5 (B) 


3 




1 


1 (D) 


1 (D) 


2.3 (B) 


3(C) 


3.7 (B) 


4.3 (B) 


5(C) 


-4 




1 


1 (D) 


2(C) 


2.7 (B) 


3.5 (B) 


4.2 (B) 


5 (C) 


5.8 (B) 


-5 




1 


1 (D) 


2.3 (B) 


3(C) 


3.8 (B) 


4.6 (B) 


5.4 (B) 


6.2 (B) 


6 




1 


1 (D) 


2.4 (B) 


3.2 (B) 


4(C) 


4.8 (B) 


5.7 (B) 


6.5 (B) 



Table 1: Numerical values on Cv in terms of units of Boltzmann's constant, k, for the 
Padmanabhan model. 



We will now treat the case with a potential of the form J2J ■ This calculation 
is rather similar to case A above. The result is 

for all D(> 2). That is, the specific heat is positive. 

We will later on be interested in the distribution of mass (or probability 
distribution) as a function of r, p r (r) (= r D ~ 1 p(r), where p(r) is the ordinary 
density or probability distribution). Observe that, for potentials of the form 0J, 
p r is the same as the integrand in J§J). The series expansion of p r (r) at r = is 

D-l 

Pr(r) = r D - x ]T hr vi (18) 

i=0 

where fe; ^ 0. For positive v, the term with the lowest power of r is an r D ~ x 
term, for v = —1 it is a constant term, and for v < —2 it is an r ("+ 1 )( D - 1 ) term. 

2.2 The Lynden-Bell model 

This model was first presented in Lynden-Bell and Lynden-Bell 5 . N parti- 
cles are, in three-dimensional space, confined to a spherical surface of radius r. 
When it comes to the gravitational interaction, the total mass of the particles, 
M, is assumed to be uniformly distributed over the surface. The radius, r, is 
fluctuating as a result of fluctuations in the distribution between potential and 
kinetic energy in the system. 

Here, we generalise the model to a space with D (D > 2) dimensions, assum- 
ing the sphere to have one dimension less than the space has. We also generalise 
the potential to be of the form 0J, where C and v have the same sign. The 
Lagrangian of the system is 

Mr 2 "^r 2 A 2 J% C „ 
< - \ \ M —r (19) 



^ ^ 2 2 

j=i »=i 



6 



where the 9ji'.s are angular velocities of particle j, the Ajf.s constants (inde- 
pendent of r, r and the 9ji'.s, but dependent of the Ojf.s). The Hamiltonian 
is 

2 N D-l 2 r , 

The phase-space volume is 

g{E)oc &r{E-M-r v )- A ^ 11 -^r N ^ D -^AE (21) 

J a 2 

where we have omitted constants (not dependent on r or E). The number r rnax 
represents the limit on r set by the total energy of the system, E. The number 
a represents an inner cutoff that keeps the phase-space volume, g(E), finite. We 
let a — > 0. We then have 

g(E) oc &EE^ N{D - 1)+ ^ v>-2 (22) 

and, again using (|TU)l 

Cv = k( H^ N ^ D - l) + v -~ 2 (23) 

We see that the specific heat is negative for v — — 1 for any D (D > 2). There 
is also a very small amount of negative specific heat, Cy = — k, for v = — 2 for 
any D (D > 2). 
We also have 

= i/ < -3 (24) 
When considering a potential of the form J5J), we get 

Cv = M ffizl)_i, (26) 

for any D (D > 2). That is, the specific heat is positive. 

We will later on be interested in the distribution of mass (or probability 
distribution) as a function of r, p r (r) (= r D ~ l p(r), where p(r) is the ordinary 
density or probability distribution). Observe that, for potentials of the form 
JTJ, Pr is the same as the integrand in (|21() . The series expansion of p r (r) at 
r = for negative v is 



Pr(r)=r^ N ^- ! iJ2 b ^ Ul ( 26 ) 



_ ^2^^(23-1)- _ 

i=0 

jy(o-i) i i 



where 7^ 0. The term with the lowest power of r is for u = — 1 an r 2 '2 
term, for v = —2 an r 1 term, and for v < —3 an r - ^ - ^^ -1 ^^ term. For 
positive ^ the corresponding expression involves only positive powers of r. 

3 Discussion 

For potentials of the form J2J), the specific heat is greater than zero, so here we 
focus on discussing potentials of the form (JTJ. As can be seen in fig. 1, there is 
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a large region where long-range forces define the system, and where the system 
has positive specific heat. Therefore, it seems to be something more than the 
long-range nature of forces that is necessary for negative specific heat. For the 
Virial theorem, the region of negative specific heat coincides with the region 
of negative total energy of the system. For the Padmanabhan model and the 
Lynden-Bell model, it is possible to chose to study negative or positive energy 
of the system. We chose to study the case where there may be negative specific 
heat, that is, the case where no outer cutoff affects the system, and the energy 
has to be chosen negative for negative v. In these cases, where all combinations 
of D and negative v are studied with negative energy, significant negative specific 
heat is still present only for v = — 1, (but not for the Padmanabhan model when 
D = 2, where the specific heat is positive). In the Padmanabhan model and the 
Lynden-Bell model, the effective potential felt by a particle due to all the other 
particles, is of the form </> = Cr", where r is the radial coordinate of the particle 
relative some fixed point, or centre of mass of the system. It would then be 
possible to suggest, that if the long-range nature of a force is sufficient to result 
in negative specific heat, this arises from emerging density distributions in the 
system that arise because of the long-range nature of the force, and which give 
an effective potential on another form than the one investigated here. These two 
models would then not be suitable for this kind of analysis. The results obtained 
from them are, hovewer, in agreement with the result from the Virial theorem. 
The Virial theorem relies on the mutual interaction between pairs of particles of 
the form Cr v , and is valid independently of any emerging density distribution. 
This then becomes an argument that there has to be something more than the 
long-range nature of the force to give negative specific heat. To regard the 
stability of a system and the Virial theorem to be this fundamental mechanism, 
does not shed sufficient light on the dynamical aspects of negative specific heat. 
Furthermore, the Virial theorem applied to systems with interaction potentials 
that are not homogeneous in r, does not give a direct relation between the kinetic 
and the potential energy of the system. The radial density (or probability) 
distribution then enters the formula. This function is usually difficult to derive. 

As mentioned before, long-range forces give rise to non-extensivity. In 
Thirring [5], a model with a potential on a quite another form than the one 
studied here is presented. In his model it is possible to vary a parameter to 
change the degree of non-extensivity of the energy of the system. Thirring's 
analysis shows that negative specific heat enters near to the point where the 
system becomes non-extensive. 

The definition of long- and short-range forces used here might not be the 
relevant one. Since the systems studied are not homogeneous, other definitions 
not relying on a matter distribution with constant density may be more accurate. 

When it comes to the very important question of a physical mechanism 
behind negative specific heat, we will here try to give an interpretation that 
is consistent with the results of the two soluble models, and that is intuitively 
satisfying. It is built on the observation that systems with particles mainly 
affected by walls, exhibit positive specific heat. With "walls", we here mean 
any potential that keeps particles in a limited volume of space, so that the 
specific heat is positive. This also includes "smooth" walls, like a potential of 
the form JQ) with v = 1 , which is sufficient to keep the particles from escaping to 
infinity, and that gives positive specific heat according to the two soluble models 
and the Virial theorem. All positive v represent this kind of wall, which we call 
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a wall of the first kind. This kind of wall may roughly be characterised with 
the criteria <f>"(r) > 0. There is a second kind of wall, roughly characterised by 
4>"{r) < and having both a very steep part (representing the wall) and a flat 
part (enabling negative specific heat). One example of this kind of wall is the 
very steep potential at small r when v is negative. To result in positive specific 
heat for a complete system, this kind of wall requires that many of the particles 
in the system are close to each other. This is clearly true for systems where the 
density or probability distribution p r (r) (= r D ~ 1 p(r), where p is the common 
density or common probability distribution) has a singularity for r — 0. (The 
maximum density will in our models occur at r — a, where a is very close to 
zero.) This requires an explanation: For not very small r, the density p is a 
measure of how many particles that keep each other close with walls of the 
second kind. But with a very sharp maximum at r = a, and a very small, all 
particles at r = a will keep each other close. Therefore p r is here the correct 
measure. Observe that p r is the same as the integrand in © and (|21J) . The 
series expansions of p r at r = is expressed in (fl%|l and (|2l))) . For those cases 
where we have negative v and positive or zero specific heat, p r has a singular 
point at r = 0. Negative specific heat only appears for those cases where v is 
negative, and there is no singular point at r = 0. The only exception from this 
principle is the Padmanabhan model in the case with D = 2 and v = — 1. For 
the Lynden-Bell model with v = — 2, the negative specific heat is very small, 
but the sign agrees with the principle. For a potential of the form (J2J, we also 
have walls of the second kind. The principle in this case holds for all D in both 
models, since the specific heat is positive and p r always has a singular point at 
r = 0. The principle described above is independent of the Virial theorem. It 
is built upon properties of the potential, and upon characteristics of emerging 
density (or probability) distributions. (Observe that the same mechanism would 
probably work even if there were sharp peaks in p r instead of singular points.) 
The principle can also be used to explain the results of Thirring |5] when the 
energy is a superextensive parameter. A potential well surrounded by a domain 
with infinite potential energy confines the particles, independent of how high 
their energies are. At low energies, the particles are confined in a smaller, even 
deeper well inside the former one. Both at high and at low energies, the system 
is then mainly affected by walls of the first kind. This agrees with the results 
of Thirring that show positive specific heat in these intervals of energy. In 
an intermediate interval of energy, the system is mainly affected by a wall of 
the second kind for which, in the spherical symmetric case, there is a region 
in r where 4>" {r) is negative and big. (We here assume an approximation to 
the potential, such that second derivatives with respect to space parameters 
are defined everywhere.) Since the potential is not steep at short distances 
between particles, it is obvious that in this system no singular points in any 
density distribution can be expected, and no other ways for the particles to 
create subsystems with many particles and positive specific heat, seem to exist. 
This also agrees with the results of Thirring, which show negative specific heat 
in this domain of energy. 
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4 Conclusions 



For systems consisting of particles interacting with potentials of the form 
the Padmanabhan model and the Lynden-Bell model give the same interval 
(in D and v) of negative specific heat as the Virial theorem does, for D > 3. 
Significant negative specific heat is present for v = —1. Long-range forces obeys 
D + v > 0. Hence, the long-range nature of a force is not sufficient to give 
negative specific heat. 

To use the Virial theorem to identify negative specific heat in systems with 
arbitrary interaction potentials, the radial density (or probability) distribution 
has to be regarded. To derive these distributions, and their relation to energy 
and temperature, is crucial to identify negative specific heat in these systems, 
if the Virial theorem will be used. 

A physical mechanism behind negative specific heat, independent of the 
Virial theorem, is suggested. It implies that positive specific heat appears when 
the system is strongly affected by what we call walls of the first kind, or the 
density distribution shows that many particles are affected by what we call walls 
of the second kind. In all other cases negative specific heat is assumed to be 
present. This is consistent with the results of the two soluble models, except for 
in one case (the Padmanabhan model with D — 2 and v = —1). The principle 
may, in some developed form, be an important part of "the general cause to 
negative specific heat" . 
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. | Gravitation 



X Negative specific heat (for both models and the Virial theorem) 

Negative specific heat (for the Lynden-Bell model and the Virial theorem) 

O Zero specific heal (for the Virial theorem) 

and approximately zero specific heat (for the Lynden-Bell model) 



f^J Zero specific heat (for the Padmanabhan model) 



All other cases involve positive specific heat 

except for the Lynden-Bell model thai involves zero specific heal for v less than or equal to -3 
A: Walls of the first kind 
B: Walls of the second kind 

C: Singular points (except for the Lynden-Bell model with v = -2) 



Fig. 1: For potentials of the form Q: Existence of negative specific heat as a func- 
tion of D and v for the Virial theorem, the Padmanabhan model and the Lynden-Bell 
model. Existence of singular points in p r for the Padmanabhan model and the Lynden- 
Bell model. Existence of walls of the first kind and of the second kind. 
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